Spin-orbit induced spin-density wave in a quantum wire 
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We present analysis of the interacting quantum wire problem in the presence of magnetic field 
and spin-orbit interaction. We show that an interesting interplay of Zeeman and spin-orbit terms, 
facilitated by the electron-electron interaction, results in the spin-density wave (SDW) state when 
the magnetic field and spin-orbit axes are orthogonal. This strongly affects charge transport through 
the wire: with SDW stabilized, single particle backscattering off an nonmagnetic impurity becomes 
irrelevant. Sensitivity of the effect to the direction of the magnetic field can be used for experimental 
verification of this proposal. 
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Introduction. Active current interest in devising 
schemes to manipulate electron spin has led to several 
interesting developments 0, 0, [j| • Most approaches rely 
on spin-orbit interaction, which couples particle's mo- 
mentum and spin, in order to achieve this goal. While 
significant progress has been made in clarifying the role 
of spin-orbit (SO) interaction, mostly of the Rashba type 
[3], on the electric and spin transport by nonintei acting 
electrons during the past few years |5(, our understand- 
ing of the combined effect of SO and electron-electron 
interactions is still limited [1, 0, B 0] • 

Here we study the combined effect of (Zeeman) mag- 
netic field and spin-orbit interaction in a single-channel 
interacting quantum wire. This set-up allows for the 
well-controlled theoretical analysis of the interplay be- 
tween broken time reversal T (by applied magnetic field) 
and inversion V (by spin-orbit interaction) symmetries 
and electron-electron interactions. The problem is for- 
mulated as follows. We consider a single-channel bal- 
listic quantum wire, corresponding to the two-terminal 
conductance Go = 2e 2 /h. Applied magnetic field creates 
two spin-split subbands, the wave functions of which are 
given by the standard orthogonal pair, (j | = (1,0) and 
(i | = (0, 1) (the orbital effect of the field is neglected). 
It reduces spin-rotational symmetry to U(l), rotations 
about o z axis. Next we add weak spin-orbit interaction 
= a RPy®x, which is obtained by electrostatic gat- 
ing of two-dimensional electron gas with Rashba SO in- 
teraction [3j|. (Corrections to this form, due to the omit- 
ted "transverse" piece, a^p x a y , and virtual transitions 
to the higher, unoccupied, subbands, can be taken into 
account [6| but are irrelevant for our purposes here.) Ob- 
serve that -ffjj 1 ^ breaks spatial inversion (y — > —y) and 



[7(1) spin symmetry, [<j z , H jj lc ^] ^ 0. The major con- 
sequence of this is the ope ning of a new, inter-subband 
Cooper scattering channel [l(| [n[. In this process, a 
pair of electrons with opposite momenta in one subband 
is scattered ("tunnels") into a similar pair in the other 
subband, see Figure [TJ Note that this process requires 
spin non-conservation (i.e. ur ^ 0), mentioned above, 
as it scatters two electrons with (almost) "up" -spins into 



a pair with (almost) "down" -spins (and vice versa). This 
simple observation is the key to our analysis: its deriva- 
tion and consequences are presented below. 

Technical formulation. The single-particle Hamilto- 
nian, describing the scenario outlined above, reads 

p 2 1 

Hq = - p - -gpL B Bu z + a R po- x , (1) 

2m 2 

where momentum along the wire (y-axis) is now denoted 
as p. The eigenstates ip u (y) = e %vv \ v (p) (v = =p) are 
easily expressed in terms of the momentum-dependent 
spinors [2. 1X2 1 
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(2) 



which describe momentum-dependent orientation of elec- 
tron's spin in z — x plane. The rotation is specified by the 
angle 7 P = arctan(2a rp / gfisB)- Note that left- (right-) 
moving particle experiences clockwise (counter-clockwise) 
rotation of its spin away from "up"-spin (v = —) and 
"down"-spin (v = +) orientations at the subband's cen- 
ter p = 0, see Figure [TJ The corresponding eigenvalues 

2 , 

Ezp = — pT v (mrp) 2 + (gn B B/2) 2 describe two non- 
intersecting branches. The gap between them is again 
momentum-dependent and is minimal at p = 0, where it 
reduces to the Zeeman energy gfisB. 




FIG. 1: Occupied subbands E^. of Eq. [TJ Arrows illustrate 
spin polarization in different subbands. Dashed (dotted) lines 
indicate exchange (direct) Cooper scattering processes. 
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We consider the situation when the Fermi energy 
Ep = vfPf, where pp (vf) are Fermi-momentum (ve- 
locity), crosses both branches, as shown in Fig. [TJ result- 
ing in four Fermi points, ip^ in the wire. To describe 
low-energy excitations of the interacting wire we project 
single-particle spin-s state ^ s (s =T,|) onto the two- 
dimensional space spanned by tp^ eigenstates: 



{Xu{p v )\sy™R v + { X u{-Pu)\s)e-^yL v . 



(3) 

Operators R v [L v ) represent slow degrees of freedom: 
right- (left-) movers in the vicinity of +p u {—pv) Fermi- 
points of the v-ih subband, respectively. In this represen- 
tation, the interaction term H mt = ^ ^ s s , J dydy'U(y — 

y') x Sl(y)w s ,(y')^! s /(y / )^! s (y) reduces to the sum of intra- 
and mter-subband scattering processes Keeping 
only low-energy momentum-conserving ones, the inter- 
subband terms include, in notations of Ref. I 111 , forward, 
exchange-backscattering and Cooper processes. The 
Cooper scattering represents two-particle (pair) tunnel- 
ing between — and + subbands. It reads 

Hc = Jdy {U(p- - P+ ) sin 2 [( 7 _ - 7+ )/2] (4) 

-U(p- + p+) sin 2 [( 7 _ + 1+ )/2]}{rIlIr + L+ + h.c.) 

Here U(q) — J drU(r)e lqr is the q-th Fourier compo- 
nent of electron interaction. The terms inside the figure 
brackets in (j4]) represent matrix elements for two differ- 
ent Cooper scatterings - direct and exchange, see Fig. Q] 
U(p- — p+) describes direct scattering in which right- 
mover R v in the z/-th subband scatters into right-mover 
R-v in the opposite, —v, subband, R u — > while 
its left-moving companion L v scatters into L_, y . The 
other possibility, exchange Cooper scattering, involving 
U(p- +p+), describes right and left members of the pair 
scattering across: R u — > L- v and L v — > R- v . It is crucial 
to observe here that in addition to involving two differ- 
ent Fourier components of the interaction potential, these 
two processes include squares of single-particle overlap 
integrals, sin 2 [( 7 _ =p 7+ )/2]. Relative magnitude of these 
is easy to understand in the limit of strong magnetic field 
and weak spin-orbit splitting, anp±/ '(g/is-B) <C 1, on 
which we concentrate now. As discussed in the introduc- 
tion, in this limit eigen-spinors Xt a l m ost coincide with 
spin \s =Tj!) eigenstates of the Zeeman Hamiltonian. 
Weak SO term, which can be thought of as momentum- 
dependent magnetic field, acting along the orthogonal, 
& x , direction, causes spins at p- and p+ Fermi-points 
to tilt by only slightly different amount, resulting in a 
small overlap of single-particle wave functions, propor- 
tional to the difference 5pF = P- — P+ = g^BB /vf- 
At the same time, spins at, say, right p- and left —p+ 
Fermi-points, tilt in opposite directions, resulting in rel- 
atively large angle (and bigger overlap) between them, 



proportional to the sum p_ + p+ = 2pp. This al- 
lows us to estimate the ratio of the two amplitudes as 
(U{5pF)/U(2p F )){gHBB/EF) 2 < 1 and neglect the con- 
tribution of the direct Cooper process in the following. 

Bosonization. We now bosonize the problem [l3| with 
the help of two conjugated fields, ip v and Q Vl obey- 
ing commutation relation [tp u (x), 9 v i (y)] = (i/2)6„„'(l — 
sign(a: — y)). Fermions are represented as R v — 
r\ v exp ^y^F ^ — 8 v )]/V2~7Ta and L v — r\ v exp[— iy/n((p v + 
Qv)\l \]2ita. Klein factors r\ v , satisfying {rj u , = 25 vu ' , 
insure anticommutation of fermions from different sub- 
bands, and a ~ p F x is a short-distance cutoff. We then 
transform to convenient symmetric, ip p — (ip- +ip + )/ \[2, 
and antisymmetric, (p a = (tp- — ip + )/y/2, combinations 
(and similarly for p / a ), in terms of which the Hamil- 
tonian of the problem decouples into two commuting 
ones. As indicated by notations, symmetric (antisym- 
metric) combinations in fact coincide with the standard 
charge (spin) ones. This is not a generic property of 
the problem but rather a convenient feature of the limit 
olrPf <C g/is-B <C Ep which is used in the rest of the 
paper. Symmetric (charge) part H p is purely harmonic 

I 



H P = 



K, 



Hdy^p) 2 + v F (d y e p f 



(5) 



with stiffness K p 1 



yj\ + (2*7(0) - U{2p F ))/irv F - The 
antisymmetric (spin) one includes nonlinear cosine term, 
representing Cooper process Q 

' r ^(dy^r+VFidye^ 
V 



H a = 



■ cos\v8tt9, 



K. 



- 1 = y/1 - U(2p F )/7rvF, 9c = U(2 PF )[ 



( 2a R p F \ ' 



g[i B B J 

Renormahzed velocities of these excitations follow from 
v p,a = vp/Kp^. Equations (|5I6|) include H (fT|) 
as well as momentum-conserving intra-subband (for- 
ward and backscattering) and inter-subband forward (cx 
U(0)) interactions, which are encoded in the stiffnesses 
K p j a . Inter-subband exchange backscattering, although 
momentum-conserving, is neglected because it is strictly 
marginal and small, of the order a R . We have also omit- 
ted marginal correction, small in gugB/Ep <C 1 fac- 
tor, associated with weak dependence of subband veloc- 
ities Vzf on magnetic field 1J] - this is the main reason 
for the equivalence of symmetric (antisymmetric) modes 
with charge (spin) ones, mentioned above. Yet another 
simplification consists in replacing U(2p±) by U(2pf) in 
expressions for K p / a - this is a valid approximation for 
any physical U(r). Finally, we must keep Cooper term 
in ©, which, inspite of having small amplitude g c , is 
strictly relevant in the renormalization group (RG) sense. 
Its scaling dimension is 2/K a < 2 for repulsive interac- 
tions jlB ]. 

Full argument in favor of Cooper term's relevancy is 
a bit more delicate. It has to do with irrelevant inter- 
subband direct backscattering term tx g^s cos[\/8tt ip^ — 



(6) 
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2Sp F y], omitted from ((6]). Note that K a = 1 + g^/2. 
Backscattering decays as ghs{£) = 5bs(0)/(l + gb s (0)£) 
until the rescaled cutoff reaches a{£) = ae £ ~ l/Sp F , see 
0. At that scale t = \n(p F /5p F ) = ln(E F /g/j, B B) 
and strongly oscillating spin backscattering cosine dis- 
appears from the problem ("averages out") [l5j . Spin 
stiffness K* = 1 + g^ s {£*)/2 stops at the value above one 
(lij ]. which implies the relevancy of the Cooper term, as 
already mentioned above. In more detail, the Cooper 
coupling constant, evolution of which is described by 
the simple dig c = (2 — 2/K a )g c , changes little from its 
initial value by the time scale £* is reached: g c (i*) = 
<7 C (0)[1 + <7bs(0)^*]. From this point on, one is allowed to 
neglect g\, s completely, and treat the Cooper scattering 
term Eq. [6] as the only relevant interaction. Both g c and 
K c grow under RG as £ is increased past £*, and reach 
strong coupling limit when g c (£) ~ v F while K a — > 2 [l6| . 

Consequences of |6p. The flow to strong coupling im- 
plies the change in the ground state (of spin sector) from 
gapless to gapped. The resulting spin gap can be esti- 
mated as A - [{a R p F / gn B B) 2 U(2p F )/v F ) K °IW K °- 1 » . 
This gap represents energy cost of (massive) fluctuations 
89 a near the semi-classical minima Q a = (jn-\-\)-\f%, m £ 
Z of #-field. Physical meaning of these minima fol- 
lows from the analysis of spin correlations. Choosing 
the gauge where r)^r]i = 
components of spin 



S x 
S z 



17| . we find for the 2p F - 



cos[V2TTip p + 2p F y] 



2pn 



sinfv^TT^o-] 
- cos[V2tt9 
sm[\/2nip cr ] 



cos[v27r^p + 2p F y] 



±1 





(7) 



The last line of the equation above is somewhat symbolic, 
with zeros representing exponentially decaying correla- 
tions of the corresponding spin components, S v,z . Here 
y-component does not order because cos[\/2~tv0 = 0, 
and z-component is disordered by strong quantum fluc- 
tuations of duai if a field, as dictated by [ip, 9] commu- 
tation relation. Thus "Cooper" order found here in 
fact represents spin-density- wave (SDW^) order at mo- 
mentum 2p F of the i-component of spin density. Ob- 
serve that S x ordering is of quasi-LRO type as it in- 
volves free charge boson, ip„. As a result, spin correla- 
tions do decay with time and distance, but very slowly 
(S x (x)S x (0)} ~ cos[2p F x] x~ K ". As K p < 1 in interact- 
ing quantum wire, this is slower than x~ 1 decay typical 
for a one-dimensional Mott insulator - antiferromagnetic 
Heisenberg chain [l3T |. 

There is one more, very intriguing, consequence of 
SDW X order: suppression of 2p F charge fluctuations. 
Indeed, 2pi?-component of the charge density operator 



cos [v 271-00-] 







(8) 



reads, keeping subleading (oc 7 PF ) contribution 

p(yh PF = sm[V2irtpp + 2p F y] ( cos[V2w ip. 

2anp F 
gfi B B 

The first term is standard, and represents intra-subband 
contribution, while the second, involving 9 ai is due to 
the subleading inter-subband contribution, which couples 
± bands. Observe that both contributions disappear in 
the SDWz phase (9 a — > 9 a ). Since 2j?F-component of 
the charge density describes backscattering (p — > —p) 
of electrons by potential impurity, Eq. [8] implies irrele- 
vancy of the impurity in the SDW state. The reason 
for this is somewhat similar to that of backscattering 
suppression in the spin Hall effect |l8|: in SDW K phase 
right- and left-movers within a given subband have op- 
posite (orthogonal) S^-components, as can be seen from 
([7]) and Fig. [TJ which forbids intra-subband backscat- 
tering. (In the spin Hall case right- and left-movers 
form Kramers pair and backscattering is forbidden by 
the T symmetry [lj|, which is broken here.) Figure [I] 
also suggests that backscattering between right-movers 
of the + subband and left-movers of the — one is possi- 
ble: their S x components are parallel. Nonetheless such 
backscattering is still suppressed because of the destruc- 
tive interference of the two scattering paths. Namely, 
the inter-subband part of 2p F density oscillation, bosonic 
form of which is given by the second term in ([5]), reads 
(R_L+ — Rill- +h.c.) in terms of original fermions. The 
crucial relative minus sign between the two backscatter- 
ing processes can be traced to Eq l2l3l and represents the 
noted destructive interference. It is useful to understand 
this result perturbatively: the intra-subband piece of (JHJ) 
arises from fusing tp a from the localized impurity poten- 
tial (first term in ([SJ) with that in H^ d \ This explains 
its magnitude (oc ct B p F /5p F ) and oddness under inver- 
sion (about the impurity site) V . Thus, a potentially 
more relevant, but even under V , backscattering process 
{R ] _L + + R) + L- + h.c. ~ cos[V27T0p] sin[\/27T6' (T ]) (note 
the relative plus sign) can not be generated. 

Although the single-particle backscattering 
pressed, the two-particle in general is not 
By considering fluctuations 59 a one indeed generates two- 
particle backscattering term oc (V 2 / A) cos{V&tt ip p ]. This 
spin-insensitive impurity affects finite-temperature linear 
conductance as G-2e 2 /h oc -{V 2 / A) 2 T iK »~ 2 [11]. The 
correction is seen to become strong (relevant) for strongly 
interacting wire with K p < g, when the impurity cuts off 
charge transport completely [2l[. This leaves us with 



the finite window, ^ < K p < 1, where the impurity is 



irrelevant. This is an interesting, and, to the best of 
our knowledge, new, addition to the Kane-Fisher result 
of always relevant impurity in a single-channel repulsive 
Luttinger liquid 2l|. Note however that our discussion 
assumes fully developed SDWa; phase, and thus implies 
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weak disorder potential V <C A. Complete solution re- 
quires simultaneous RG analysis of Cooper and impurity 
terms (l7j . 

The correlated state can also be probed via tunneling 
density of states (DOS) measurements. Skipping details, 



which are rather similar to the calculation of DOS in [ll| . 
we quote the result for the local DOS in SDW X state: 
u(w) oc Q(w - A)(w - A) fc , where b = (K p - l) 2 /{4K p ) 
and O denotes the step function. Naturally, DOS is zero 
for energies below the SDW gap, and is found to rise 
smoothly (b > 0) just above it. 

Angular stability of the SDW X state can be ana- 
lyzed via angular dependence of subband dispersions 
E T in Fig. [TJ Indeed, suppose that the two axes, B 
and SO, are not orthogonal and denote the angle be- 
tween them as tt/2 — (3. This will modify the SO 
term in (fTJ) to aup(a x cos (3 + o~ z sin /?) . The eigenval- 
ues of the modified Hamiltonian (fT]) now read E T = 

- M T J K cos f3p) 2 + {\g^ B B - a R sin (3p) 2 , and 
describe two subbands (=f) shifted in the opposite di- 
rections along the momentum axis. For small [3 the dis- 
persion can be approximated as E T = ^ P %^ M T 

\J(otRp) 2 + (gpsB/i) 2 . Thus the lower (— ) subband 
shifts left and is centered around —po, while the up- 
per (+) one shifts toward positive momenta and cen- 
ters around +po, where po ~ man(3. This simple ob- 
servation implies that opposite-Fermi-momenta pairs in 
± subbands acquire opposite (±Po) center-of-mass mo- 
menta. This can be pictured by shifting the bands in 
Fig. [1] horizontally in opposite directions. Thus two- 
particle Cooper tunneling processes illustrated in Fig. [T] 
become momentum-non-conserving ones. As a result, 
this important scattering channel will disappear above 
some critical misalignment angle j3* which can be esti- 
mated as follows. Cooper order is destroyed once the 
misalignment cost oc 2vfPo becomes comparable to the 
Cooper gap A. Estimating the latter at K a — 2 we find: 
/3* w a R p F U(2p F )/(gp B B) 2 <C 1. This estimate shows 
that the found SDW^ has narrow but finite region of 
angular stability, and agrees fully with results of more 
detailed RG-based calculations in [13]. SDW X state can 
also be destroyed by reducing magnetic field strength be- 
low the critical, gfiF)B c ~ apipF, even while maintaining 
the orthogonal orientation (angle (3 = 0). This happens 
due to the decrease of the spin stiffness K a below 1 (so 
that the scaling dimension of Cooper term © exceeds 2) 
once the Zeeman energy becomes smaller than the spin- 
orbit one [l3| ■ This weak-field region, which includes the 
B = limit of fl}, has been studied previously @,0| and 
contains no relevant Cooper processes. 

The sensitivity of the described SDW^ phase to the 
mutual orientation and magnitude of the magnetic and 
SO terms can be exploited in experimental searches of the 
novel field- induced SDW phase of the quantum wire with 
spin-orbit interaction. It appears that lateral quantum 



wells at vicinal surface of gold, which possess spin-orbit- 
split and highly one-dimensional subbands, [22j can serve 
as a nice experimental starting point. 
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